Scanning Tunneling Microscopy and Spectroscopy study of charge inhomogeneities in 

bilayer Graphene 
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We report room temperature scanning tunneling microscopy and spectroscopy study of bilayer 
graphene prepared by mechanical exfoliation on Si02/Si surface and electrically contacted with gold 
pads using a mechanical mask. The bulk conductivity shows contribution from regions of varying 
electron density indicating significant charge inhomogeneity Large scale topographic images show 
ripple like structures with a roughness of ~lnm while the small scale atomic resolution images show 
graphite-like triangular lattice. Local dl/dV-V tunnel spectra have an asymmetric V-shape with the 
minima location showing significant spatial variation indicating inhomogeneity in electron density 
of order 10^^ cm— 2. The minima in spectra at a fixed location also shifts linearly with the gate 
voltage with a slope consistent with the field induced carrier density. 



Graphene has been a focus of research for past several 
years for its fascinating physics [1[ and application poten- 
tial arising from its unusual mechanical, optical, chem- 
ical, and electrical properties. The role of disorder and 
local electronic properties on bulk transport properties 
and the Dirac nature are some of the outstanding issues. 
The disorder is inevitable as one cannot have long range 
order in 2D and also due to other effects such as local 
Coulomb impurities in the supporting substrate and the 
electron-phonon interaction. 

Scanning tunneling microscopy and spectroscopy 
(STM/S) is the most suitable tool for investigating the 
electronic nature of such defects as it can probe the lo- 
cal electronic density of states (DOS) by tunneling spec- 
troscopy by measuring local dl/dV-V spectrum [2]. Re- 
cent attempts by some research groups on graphene have 
shown a variety of results. Crommie's group has been 
able to align the standard four probe Graphene Hall bar, 
made using the electron lithography techniques, with the 
STM tip. In these samples the atomic resolution im- 
ages have been observed; however, the tunneling spectra 
do not agree with the Dirac behavior and show an en- 
ergy gap of ^lOOmeV [3] at Fermi energy. STM/S work 
on thermally grown graphene and few layer graphene on 
SiC [J-lSl also shows an energy gap of similar magni- 
tude. Such graphene is known to be strongly coupled 
to the underlying SiC substrate. Another similar study 
by Deshpande et. al. does not show such a gapped be- 
havior [7[. Andrei's [8] group was able to find a small 
graphene on bulk graphite to show all the essential fea- 
tures of graphene including Dirac nature and Landau lev- 
els but without any control on electron density. Landau 
levels have also been seen on the bulk graphite surface [9| . 
On bilayer graphene, recent STM/S work by Deshpande 
et. al. [10| shows charge impurities with no correlation 
with topography. 

In this paper, we present the STM study of bilayer 
graphene prepared by mechanical exfoliation on an insu- 
lating Si02 layer on a degenerately doped Si substrate 
acting as gate electrode. The small scale STM topo- 



FIG. 1: Left picture shows the optical image of the graphene 
flake on Si02/Si identified as bilayer using the Raman Spec- 
tra, taken with a laser light of 514nm wavelength, as shown 
in the right plot. The 2D peak in this spectra near 2700cm~^ 
can be fitted with four Lorenzians as shown in the inset as- 
certaining the flake as a bilayer graphene. 



graphic images show atomically resolved surface while 
large scale images show ripple like structures. The local 
tunnel spectra (dl/dV-V) have an asymmetric parabolic 
shape with a finite zero bias conductance and with a 
minimum whose location along bias voltage axis varies 
spatially as well as with the gate voltage. The shift in 
minimum with the gate voltage is quantitatively consis- 
tent with the field induced carrier density but with an 
attenuation factor corresponding to a screening length 
of ~1.2 nm. Both bulk conductivity and spatial varia- 
tion in local tunnel spectra are compatible with charge 
inhomogeneities. The local conductance maps have a 
week correlation with the topography implying part of 
the topographic contrast is coming from charge inhomo- 
geneities. We also find more charge inhomogeneity over 
longer length scales. 

The graphene flakes were prepared by mechanical ex- 
foliation of graphite flakes [11] on 300 nm thick Si02 
on highly doped (p <0.005r^-cm) silicon substrate. Bi- 
layer graphene sheet, size (~ 42/im x 16/im) was opti- 
cally identified [12] as shown in Fig(T^. The central large 
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FIG. 2: Gate voltage dependent conductivity of bilayer 
graphene. The black line shows the experimental data while 
the red one is the model conductivity using Eq. [21 as de- 
scribed in the text, with the parameters listed in the plot. 



portion of graphene was masked using a tungsten wire of 
25/im diameter. A thin metal sheet with a sHt of 1mm 
width and few mm length was also placed on top with slit 
length perpendicular to the wire to confine each contact 
to a narrow region. lOOnm gold film followed by lOnm 
chromium was deposited after this masking to obtain two 
electrical contacts on either side of bilayer graphene as 
shown in Figl2] inset. The graphene flake was confirmed 
to be a bilayer using Raman spectroscopy done in am- 
bient conditions after the STM/S measurements. The 
Raman spectra, taken with 2.7mW Laser light of 514nm 
wavelength, are shown in Figdl We observed the G peak 
around 1580 cm~^ and 2D band at 2700 cm~^. Inset in 
Fig d] shows zoomed-in plot of 2D band with fitted to four 
Lorentzians at peak frequencies of 2657, 2689, 2709, and 
2722 cm~^. Four peaks in the 2D band of the Raman 
S2ectrum is the characteristic feature of bilayer graphene 
[13]. The absence of D peak at 1350 cm~^ demonstrates 
the high quality of our samples. 

All STM and transport data reported here were taken 
at room temperature and in 10""^ mbar vacuum. This 
vacuum was obtained by cryo-pumping. Two probe 
transport measurements were done using a sinusoidal 
voltage of 597 Hz frequency and 100 mV amplitude across 
the graphene device in series with a 2.7 1^ resistance. The 
voltage across this resistance was measured using a lock- 
in amplifier for measuring the current. Gate voltage was 
applied to highly doped silicon electrode through a 270 
kQ series resistance. For STM measurements, one of the 
two electrodes (see inset of Figl2]) of graphene device was 
connected to the STM bias circuitry, so that the tun- 
neling current diffused in-plane through the gold film. 
For aligning graphene with the STM tip, we have used 
a homemade 2D nano positioner [14] with a homemade 




FIG. 3: Top left image shows a 690x690nm^ STM topograph 
of bilayer graphene taken at 0.4V bias and O.lnA tunnel cur- 
rent with the topographic profile along the marked line shown 
in to top right plot showing the ripples. The bottom left im- 
age (5.4x5.4nm^) shows the atomic resolution image while the 
bottom right image (35x35nm^) shows the atomic resolution 
together with large scale ripples. 



STM with a design described elsewhere [15|. This STM 
uses a commercial electronics and control software |16J . 
The electrodes defined in the previous step also served as 
guides for aligning the graphene with STM tip under an 
optical microscope. Electrochemically etched tungsten 
wire of 0.25 mm diameter were used as the STM tips. 
We also treat the etched tips with hydrofluoric acid to re- 
move the oxide layer [175. The tunneling spectra (dl/dV- 
V) were acquired using the ac modulation technique with 
modulation amplitude of 20 mV. We kept the junction 
resistance value as AGft all the local spectra taken at 
different gate voltage to keep the tip sample separation 
same for all the spectra and images. The simultaneous 
STM and STS images are obtained in constant current 
mode with the z-feedback time constant much larger than 
the time period of ac-modulation. In this mode opposite 
contrasts are observed in the simultaneous topographic 
and conductance images if the the surface is topograph- 
ically fiat but electronically inhomogeneous [18]. This 
also assumes the spectra to be monotonically rising with 
bias voltage. In our STM the bias voltage is applied to 
the sample while the tip is kept at (virtual) ground po- 
tential; thus the conductance at positive (negative) bias 
corresponds to electron DOS above (below) the Fermi 
energy. 

FigO shows the variation of electrical conductivity of 
this bilayer graphene at room temperature with gate volt- 
age. There is a prominent dip in the conductivity at posi- 
tive gate voltage of 22.4 V showing the hole doped nature 



of this bilayer graphene. The conductivity value at this 
dip is 0.37ml]~^, which is nearly twice of the minimum 
conductivity reported in graphene [1|, indicating good 
electrical contacts. Conductivity increases quadratically 
around this dip with gate voltage, until another smaller 
dip appears at a gate voltage of about -20 V. This we be- 
lieve is due to charge inhomogeneities in this large bilayer 
graphene flake. The electron density, n in the field effect 
geometry varies linearly with the gate voltage. Including 
the doping due to impurities or defects we can write. 



ed 



(Vg - Vd) 



(1) 



Here d (= 300nm) and k, {= 4) are the thickness and di- 
electric constant of Si02, respectively, and Vd is a voltage 
equivalent of the residual doping, e is the magnitude of 
the electron's charge. For no do ping , the density of states 
of bilayer graphene is quadratic [19| with energy near Ep. 
Also there is a large DOS at Ep so Ep changes almost lin- 
early with doping. Thus the DOS increases quadratically 
with doping giving a quadratic variation in conductivity 
with doping. Thus for fixed residual doping at a given lo- 



cation, the conductivity, a oc [1 + a( 
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Assuming 



a simple configuration of N different regions, with differ- 
ent residual doping, connected in series between the two 
contacts, the expression for the average conductivity as 
a function of gate voltage is given by. 
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Here ai decides the relative fraction of the region with 
doping Vi. We can qualitatively produce the features of 
the measured conductivity by using four different dop- 
ing regions with parameters as listed in Fig O Although 
this simple model captures essential features of the ex- 
perimental conductivity but it's an oversimplification and 
one needs to consider a resistor network where parallel 
regions also exist. More complications can arise if there is 
significant contribution from scattering at the boundary 
between regions. Although such charge inhomogeneities 
are confirmed by our STS spectra as discussed further, 
we cannot get the doping variation on the whole sample 
to calculate the detailed conductivity behavior. 

FiglSt^ shows the topographic images of various sizes 
with the largest area (690x690nm^) image showing rip- 
ple like structures on a lateral length scale of ~10-20nm 
and l-2nm height. The smallest scale (5.4x5.4nm^) im- 
age shows atomically resolved surface while the interme- 
diate size (35x35nm^) image shows simultaneously re- 
solved ripples and the atomic features. Simultaneously 
taken topographic and conductance images at zero gate 
voltage are shown in Fig. |4] together with some local 
spectra at the spots marked in the conductance image. 
We find that the spectra in the dark regions of the con- 
ductance image have a larger shift towards negative bias. 



i.e. Ep is above the DOS minima or in other words it cor- 
responds to more electron doping. In fact, most spectra 
at zero gate voltage show electron doped behavior with 
a shift as large as -lOOmeV in some of the spectra taken 
at several \im distance from this area. We estimate the 
spread in this shift to be ~40meV from the average shift 
of ~-40meV. This spread, from the discussion given later, 
is equivalent to the field effect doping by a gate voltage 
of '^35V, which from Eq. [T] corresponds to an electron 
density variation of ~10^^ cm~^. 

Other than this shift in minima the overall nature of 
the spectra remains same and thus we conclude that 
the conductance image contrast is predominantly arising 
from the variation in local electron density. Since these 
spectra are taken at room temperature, we do not see 
any gap due to large thermal smearing, which has been 
seen at low temperatures [10]. This kind of electronic in- 
homogeneity can be attributed to either the charge fluc- 
tuations due to substrate [20] or the ripples [7]. We also 
find some correlation between the conductance and topo- 
graphic images at small scales ('^lOOnm) which makes us 
believe that some of the contrast in the topographic im- 
age is coming from this doping variation. We have found 
a cross-correlation coefficient in simultaneous STM and 
STS images as large as -0.2 and for a few scans we have 
also seen a positive cross-correlation coefficient. This in- 
dicates that there are generic topographic features, i.e. 
ripples, which could also be contributing to the electronic 
contrast. 

Fig. [5^ shows the gate voltage dependence of the local 
dl/dV-V spectra with dl/dV being proportional to the 
DOS at energy eV with some thermal smearing. These 
spectra were acquired at the same point of an image (not 
shown). The spectra show the same shape except for a 
shift in the minima. The exact location of the minima 
was obtained by fitting a quadratic function in a small 
bias range (±100mV) about the minima as shown in Fig. 
[Sb insets. This gives us the location of Ep with respect 
to the minima as shown in Fig. [Sb insets. The variation 
of the Ep with gate voltage is shown in FiglSb; which 
shows a linear variation with gate voltage. The linear 
fitting gives a slope of (1.24±0.1)xl0~^. 

To calculate this slope theoretically for a bilayer 
graphene with 2-fold spin and 2-fold valley degeneracy, 
the Fermi wave-vector is given by, kp = v^mr, with 
n as the electron density. Assuming a tight binding 
parabolic band [1] near E^?, i.e., ek± = h Vpk'^/t±j we 
get Ep = ±fi VpTiTr/t^^ with t± (= 0.4eV) as the inter- 
layer coupling parameter and vp (=lxlO~^ m/s) as the 
Fermi velocity in single layer graphene. Using the field 
induced carrier density from Eq. [T]we get. 



From this, Ep {in eV)=a{Vg — Vd) with a 
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FIG. 4: The topographic (top left) and conductance (top 
right) images (121xl21nm^) of bilayer graphene taken simul- 
taneously in constant current mode at 0.4V bias and O.lnA 
tunnel current. The bottom plot shows the spectra taken at 
four different locations as marked in the conductance image. 
The inset in the lower plot shows the region close to minima 
to show the shift in the spectra. The lower spectra (with - 
0.05 offset) show a selection of spectra taken at various spots 
spanning several /xm area (to be added). 



which works out to be 2.53x10"^. The experimental 
value of the slope (see Fig. [Sb) is smaller than this theo- 
retical value by a factor of ~0.5. This, we believe, is due 
to the screening effects, which will give rise to the field 
attenuation by a factor, /3 = exp(— t^/A) with A as the 
screening length and tg as the graphene thickness. Using 
A (=1.2 nm) |21j and t^= 0.7 nm for bilayer graphene, 
we get /3=0.56. Thus the slope should be /3a=1.41xlO~^, 
which is in good agreement with the experimental value. 

In conclusion, we have done a STM/S study with 
atomic resolution of gated bilayer graphene prepared by 
mechanical exfoliation. The STS spectra have an asym- 
metric V-shape with non-zero DOS near E^? and with 
significant spatial variation in the location of minima. 
From this we find an inhomogeneity in electron density 
~10^^ cm~^. The behavior of spectra with gate voltage is 
in quantitative agreement with the field induced doping 
together with the screening effects. 
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FIG. 5: (a) Spectra at different gate voltages (V^) taken at 
the same position of an image (not shown), (b) Variation of 
Ei? as a function of the gate voltage. The two insets show a 
small portion of the spectra, taken at V^=+60 an -GOV, with 
the parabolic fitting to find the minima more accurately. The 
arrow at the bottom in the insets indicates the location of Ef 
with respect to the spectral minima. 



ACKNOWLEDGEMENTS 

We thank Mandar Deshmukh for discussions and also 
for providing us with graphite crystals. We also thank 
Rajeev Gupta and Shyam Lai Gupta for Raman Spec- 
troscopy. SKC acknowledges financial support from the 
UGC of the Government of India. Financial support 
from MHRD and DST of the Government of India is also 
thankfully acknowledged. 



[1] A. H. Castro Neto, F. Guinea, N. M. R. Peres, K. S. 
Novoselov, and A. K. Geim, Rev. Mod. Phys. 81, 109 
(2009). 



[2] "Introduction to Scanning Tunneling Microscopy", 2^^ 

ed., by C. J. Chen, Oxford Univ. Press (2008). 
[3] Y. Zhang, V. W. Brar, F. Wang, C. Giriti, Y. Yayon, M. 

Panlasigui, A. Zettl, and M. F. Crommie, Nature Physics 

4, 627 (2008). 
[4] V. W. Brar, Y. Zhang, Y. Yayon, T. Ohta, J. L. Mc- 

Chesneym, A. Bostwick, E. Rotenberg, K. Horn, M. F. 

Crommie, App. Phys. Lett. 91, 122102 (2007). 
[5] F. Hiebel, P. Mallet, L. Magaud, and J.-Y. Veuillen, PRB 

80, 235429 (2009). 
[6] P. Lauffer, K. V. Emtsev, R. Graupner, Th. Seyller, L. 

Ley, S. A. Reshanov and H. B. Weber, PRB 77, 155426 

(2008). 
[7] A. Deshpande, W. Bao, F. Miao, C. N. Lau, and B. J. 

LeRoy, Phys. Rev. B 79, 205411 (2009). 
[8] G. Li, A. Luican and E. Y. Andrei, Phys. Rev. Lett., 

102, 176804 (2009). 
[9] T. Matsui, H. Kambara, Y. Niimi, K. Tagami, M. 

Tsukada, and H. Fukuyama, Phys. Rev. Lett. 94, 226403 

(2005). 
[10] A. Deshpande, W. Bao, Z. Zhao, C. N. Lau, and B. J. 

LeRoy, Appl. Phys. Lett. 95, 143502 (2009). 
[11] www.graphite.de 



[12] P. Blake, E. W. Hill, A. H. Castro Neto, K. S. Novodelov, 

D. Jiang, R. Yang, T. J. Booth and A. K. Geim, Appl. 

Phys. Lett. 91, 063124 (2007). 
[13] L.M. Malarda, M.A. Pimentaa, G. Dresselhaus and M.S. 

Dresselhaus, Phys. Rep. 473, 51 (2009). 
[14] A. K. Gupta, R. S. Sinha, and R. K. Singh, Rev. Sci. 

lustrum. 79, 063701 (2008). 
[15] A. K. Gupta and K. -W. Ng, Rev. Sci. lustrum. 72, 3552 

(2001). 
[16] SPMIOO from RHK Technoogy, Troy, MI, USA. 
[17] L. A. Hackett and S. E. Creagera, Rev. Sci. lustrum. 64, 

263 (1993). 
[18] Ch. Renner and H. M. Ronnow, chapter H.9 in Scanning 

Probe Microscopy: Electromechanical Phenomena at the 

Nanoscale, Springer Series (2006). 
[19] M. Koshino, New. J. Phys. 11, 095010 (2009). 
[20] Y. Zhang, V. W. Brar, C. Girit, A. Zettl and M. F. Crom- 
mie, Nature Physics 5, 722 (2009). 
[21] H. Miyazaki, S. Odaka, T. Sato, S. Tanaka, H. Goto, A. 

Kanda, K. Tsukagoshi, Y. Ootuka and Y. Aoyagi, Appl. 

Phys. Exp. 1, 034007 (2009). 



